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Data samples corresponding to the isospin-violating decay D^'^ — > D'^-k^ and the decays D*^ 
Dtj, D*° D°Tv° and 0*° D°7 are reconstructed using 90.4 fb of data recorded by the 
BABAR detector at the PEP-II asymmetric-energy e"'"e~ colUder. The following branching ratios are 
extracted: r(D*+ D+7r°)/r(D*+ ^ D+7) = 0.062 ± 0.005 (stat.) ± 0.006 (syst.) and r{D*° 
7?°7r°)/r(_D*" D^^f) = 1.74 ± 0.02 (stat.) ± 0.13 (syst.). Both measurements represent significant 
improvements over present world averages. 

PACS numbers: 13.25.Ft, 13.40.Hq, 12.39.Fe 



The decay of any higher-mass cs meson into Dfn'^ 1] 
violates isospin conservation, thus guaranteeing a small 
partial width. The amount of suppression is a matter of 
large theoretical uncertainty according to most models 
of charm-meson radiative decay 0. One such model Q 
suggests that the decay D*^ Dfir^ may proceed via 
n'^-r] mixing. Even including such considerations, the ra- 
diative decay _D*+ ^ D^-f is still expected to dominate. 
The existence of isospin-violating decay modes such as 
D*~^ Dfir'-' is particularly relevant given the recent ob- 
servations of two narrow new Df meson states 0, Q • In 
particular, in contrast to the D*^ meson, there is no ex- 
perimental evidence for the electromagnetic decay of the 
Dsj{2317)^; current measurements place the branching 
ratio to D*^^ at less than 18% at 90% confidence level 
(CL) @. 

Besides the Dfir'^ and I?^7 final states, no other de- 
cay modes of the D*^ have been observed and none 
are expected to occur at a significant level. Only one 
previous observation of the decay D*^ — > Dfn^ is 
recorded in the literature, yielding a value of 0.062 
-aois (stat.) ± 0.022 (syst.) for the branching ratio 
r{D;+ DfTT^)/r{Dl+ D+7) 0. The analysis pre- 
sented here confirms this observation and provides a more 
precise measurement of this branching ratio. 

The decay D*° D°Tr°, in contrast to D*+ D+Tr°, 
does not violate isospin conservation and the world aver- 
age for the branching ratio is r(D*° D°7r°)/r(£)*° ^ 
D°-/) = 1.625 ± 0.20 0. As for the D*+ meson, the tt" 
and 7 decay modes are expected to saturate the decay 
width of the 1?*° meson. 

The results presented here are based on data recorded 
by the BABAR detector at the PEP-II asymmetric-energy 
e~^e~ storage rings. The data sample, corresponding to 
an integrated luminosity of 90.4 fb~^, was recorded at 
and approximately 40 MeV below the T{AS) resonance. 
Due to the unequal beam energies, the e^e~ center-of- 
mass system is boosted relative to the laboratory frame 
with w 0.55. The BABAR detector and trigger are 
described in detail elsewhere 0- Charged particles are 
detected and their momenta measured by a silicon vertex 
tracker (SVT) consisting of five layers of double-sided sili- 
con strip sensors and a cylindrical 40-layer drift chamber 
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(DCH), both operating within a 1.5 T solenoidal mag- 
netic field. Charged particle identification is provided 
by energy loss measurements in the SVT and DCH and 
by Cherenkov light detected in an internally reflecting 
ring imaging detector (DIRC). Photons are identified and 
their energies measured by an electromagnetic calorime- 
ter (EMC) composed of 6580 CsI(Tl) crystals. 

In the following paragraphs, the -D*^ measurement is 
described first. The D*'^ analysis, which uses similar pro- 
cedures for signal extraction, is described afterward in 
less detail. 

mesons are reconstructed via the decay sequence 
£)+ (/'TT''", (j) ~^ K^K~ . Kaons are identified by com- 
bining the energy deposited in the SVT and DCH with 
the information from the DIRC. Tracks not identified as 
kaons according to the particle identification criteria are 
considered to be pions. All K'^K~'n^ candidates are re- 
quired to fit successfully to a common vertex. Only com- 
binations with a K^K~ invariant mass within 8 MeV/c^ 
of the nominal (j) mass are retained. 

In e^e^ annihilation to charm quarks, the cc fragmen- 
tation process is characterized by the production of high- 
momentum (leading) charm hadrons. This property is 
exploited in order to reduce substantially the combinato- 
rial background by retaining only those <j)i:'^ candidates 
with scaled momentum Xp greater than 0.6, where Xp is 
defined as Xp{D+) = p*{D+)/p*^^^{D+) and p*{D+) is 
the momentum of the Df candidates in the e~^e center- 

of-mass frame with p*^^^iD+) = \J E^^^J - m{Dtf as 
its maximum value. 

The longitudinal polarization of the (j) nieson in the 
Z?j" rest frame is used to reduce background by requiring 
that the absolute value of the cosine of the helicity angle, 
defined as the angle between the cj) momentum direction 
in the D'^ rest frame and the momentum direction of 
either of the kaons in the rest frame, is 0.3 or greater. 

The resulting K^tt^ invariant mass distribution is 
shown in Fig. ^ This distribution can be modeled by 
the sum of two Gaussian functions (to represent the sig- 
nal) and a third-order polynomial (to represent the back- 
ground). The resulting binned fit yields 73 500 ± 300 
events (statistical errors only). A 13+ candidate is re- 
tained if its invariant mass is within 12 MeV/c^ of the 
nominal mass |^. 

A 7r° candidate is reconstructed by combining two pho- 
ton candidates that fulfill the following requirements. 
Each photon candidate is identified by a calorimeter clus- 
ter that is not associated with a charged track and has an 
energy in the laboratory frame of at least 45 MeV. Ad- 
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xlO^ xlO-' 




m(KKji)[GeV/c^] Am [GeV/c'] Am [GeV/c"] 



FIG. 1: The K~-r^ mass distribution. The dots represent 
data points with error bars corresponding to statistical uncer- 
tainties (these uncertainties are smaU enough that the error 
bars are difficuh to distinguish). The sohd curve shows the 
fitted function. The dashed curve indicates the background. 

candidates are defined by the region between the vertical 
dotted lines. 



ditionally, to help remove the background from hadronic 
showers, the fractional lateral width ^3|> which describes 
the shape of the shower in the calorimeter, is required to 
be less than 0.55. The fiducial acceptance of photon can- 
didates is restricted by the angular range of the EMC 
(—0.92 < cosO < 0.89, where 9 is the polar angle in the 
center-of-mass frame jyf). 

A vr" candidate is retained if it has a momentum p* in 



the 



center-of-mass frame greater than 150 MeV/c. 



Furthermore, the absolute value of the cosine of the de- 
cay angle, 6**, which is defined as the angle between the 
direction of one of the photons in the 7r° rest frame and 
the direction of the 7r° candidate in the center-of-mass 
frame, is required to be less than 0.85. For tt*' 77 
decay, the cos 6* distribution is uniform, while it peaks 
near ±1 for random 77 combinations. 

Only 77 pairs within a specified mass interval are re- 
tained. This interval is defined by the values of mass at 
which the signal portion of a function fitted to the 
77 mass distribution falls below 0.2 times its maximum 
value. This requirement accommodates the asymmetric 
shape of the 77 mass distribution and takes into account 
variations in detector calibration. A kinematic fit is ap- 
plied to the surviving 77 pairs to constrain their mass to 
the nominal tt*^ mass. 

After combining the 13+ and 7r° candidates in a 
search for the decay D*+ -D+7r°, a fit is apphed to 
the distribution of the mass difference /^m{D'^ iv^) = 
m{K+K-Tr+-K°) - m{K+K--K+). The fit function is the 
sum of a double Gaussian function to represent the signal 
and the function 



/i(Am) = iV 1 -exp 



Am — min'^) 



(Ato^ -I- aAm + b) , 



(1) 



where m(7r*') is the tt^ mass, and N, /i, a, and b are free fit 
parameters to describe the background. The exponential 
term models the kinematic threshold; this threshold term 



FIG. 2: The D*+ signals: (a) m{K+ R- tt+tt^) - 
m{K+K-TT+); (b) m{K+ R- ir+'y) -m{K+ K-n+). The dots 
represent data points. The solid curve shows the fitted func- 
tion. The dashed curve indicates the portion of the fit asso- 
ciated with background. 



has little influence on the background shape near the 
signal region. The result of this fit is shown in Fig.|2Ia). 
A signal event yield of 560 ± 40 (statistical error only) is 
obtained. 

For the reconstruction of the decay -D+7, a 

calorimeter cluster that is not associated with a charged 
track is considered a photon candidate if it fulfills the 
following requirements: the energy must be 50 MeV or 
greater in the laboratory frame and 100 MeV or greater 



the 



center-of-mass frame and the fractional lat- 



eral width must be less than 0.8. To reduce the back- 
ground due to photons from 7r° decay, a photon candi- 
date is discarded if it forms a tt'^ candidate with any 
other photon candidate in the same event. In this case, 
a 77 combination is considered a tt" candidate if the in- 
variant mass is in the range 115 < 771(77) < 155 MeV/c^ 
and if the total energy is at least 200 MeV in the e^e~ 
center-of-mass frame. 

To obtain the D*^ Df-f signal event yield, a 
fit is applied to the distribution of the mass difference 
Am(i:'+7) = m{K+K--K+-/) - m{K+K-TT+). The fit 
function is a sum of a third-order polynomial to model 
the background plus a function first introduced by the 
Crystal Ball collaboration for the signal 



/2(Am) - iV- 



exp 



(Am-/j)^ 



if (Am — j-i)/(J < a 
if (Am — /i)/cr > a 
(2) 

where iV, /i, cr, n, and a are free fit parameters and A 
and B are chosen such that the function and its first 
derivative are continuous at (Am — ^)/a — a. The fit 
result is shown in Fig.|2Ib). A signal yield of 15 600 ±200 
events (statistical error only) is obtained. 

The reconstruction efficiencies are determined using a 
Monte Carlo simulation based on 30 000 events for each 
decay mode. The simulated events are analyzed us- 
ing the same procedure as for real data. By calculat- 
ing the ratio of the number of reconstructed to gener- 
ated events, efficiencies of e(r'+7r") = 0.041 ± 0.002 and 
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TABLE I: A summary of the relative systematic uncertainties 
in the branching ratio measurements. 



.2 0.08F 



Sources 


Relative Uncertainty (%) 


V{D*+ ^ 73+^0) 






Y{D*° ^ D°^) 


Background shape 


4.8 


0.1 


Monte Carlo statistics 


5.0 


5.4 


Signal model 


3.6 


3.8 


p* dependence 


6.8 


2.8 


Quadrature Sum 


10.2 


7.2 



e{D+-i) = 0.071 ± 0.002 are found for the two D*+ de- 
cay modes. The efficiency ratio is e{DfTT'~')/e{Dfj) = 
0.58 ± 0.03 (statistical error only). 

Various sources of systematic uncertainties are studied. 
To verify that the Monte Carlo events model the data 
correctly, r decays with one or two 7r° mesons in the 
final state are studied to obtain energy-dependent Monte 
Carlo efficiency corrections for 7r° mesons and photons. 
Although this procedure indicates that no correction is 
necessary, the errors on the correction functions represent 
uncertainties in the Monte Carlo model and contribute a 
systematic uncertainty of 3.6%. 

To test for uncertainties in the background shape of 
the mass difference distributions, upper and lower side- 
bands in the K~^K~7t'^ and 77 mass distributions are 
considered. Positive signal yields are expected in these 
sidebands from either mis-reconstructed or unassociated 
7r° candidates. To measure these yields, the same fit 
functions used to determine the signal yields are applied 
to the mass difference distributions of the sideband sam- 
ples. Any discrepancy in yield so obtained from data and 
Monte Carlo simulation is considered a systematic uncer- 
tainty (4.8%). Most of this uncertainty is attributed to 
the relatively large background in the Dfi:^ decay mode. 

The measurement of T{D*+ D+7r°)/r(D*+ 
Df-f) is repeated for the subsamples of candidates within 
various p* intervals. By fitting either a constant function 
or a first-order polynomial to the branching ratio as a 
function of p* , it is possible to verify that the measured 
branching ratios are independent of p* (see Fig.El). Nev- 
ertheless, it is assumed conservatively that any p* de- 
pendence arises from unknown momentum dependencies 
of the efficiencies that may not cancel in the branching 
ratios. The difference (6.8%) between the branching ra- 
tio represented by the constant function and the integral 
of the first-order polynomial is therefore reported as a 
systematic uncertainty. 

The systematic uncertainties are summarized in Ta- 
ble ^ Combining all contributions in quadrature, a total 
systematic uncertainty of 10.2% is derived for the mea- 
surement of r{D*+ DfT:°)/T{D*+ D+7). 

The ratio r(D*o ^ D°TT")/r{D*° -> D%), where 
Z?" — > X^TT"*", is measured using the same selection crite- 
ria for the tt'^ and photon candidates as in the reconstruc- 




1.9 



1.7 



(b) 



3.8 4.2 4.6 5 
p* [GeV/c] 



3.4 3.8 4.2 4.6 5 
p* [GeV/c] 



FIG. 3: The measured values of (a) r(L>*+ 
Dtn°)/r{Dt+ Df-y) and (b) r(D*° ^ D°n°) /r{D*'' ^ 
D'^'y) in intervals of center-of-mass momentum p*. The error 
bars indicate the associated statistical error. The solid line 
(dashed line) is the result of a fit to a constant (first-order 
polynomial). The dotted line is the result from the entire 
sample integrated over p* . 



xlO' 



xlO" 




0.15 0.16 
Am [GeV/c 



FIG. 4: The D*° signals 
(b) m{K~-K'^'y) —m{K~-K 
The solid curve shows the fitted function. The dashed curve 
indicates the portion of the fit associated with background. 



(a) m{K 7r+7r°) - m{K 7r+) and 
'^'^ The dots represent data points. 



tion of D*+ D+7T° and D*+ -> D+7. To be included 
in the D'-' K^tt^ sample, a candidate and 7r+ 
combination must yield an acceptable fit to a common 
vertex and the scaled momentum Xp of the resulting 
candidate must be 0.6 or greater. Fitting the sum of a 
double Gaussian function and a third-order polynomial 
to the resulting K~'it'^ invariant mass distribution (not 
shown) produces (996.0 ± 1.5) x 10"^ signal events (sta- 
tistical error only). A K^ir^ combination is retained if 
its mass differs by less than 17 MeV/c^ from the nominal 
D" mass 8]. 

The D'^ candidates are combined with all 7r° can- 
didates; the resulting mass difference Am(Z3°7r'') = 
m{K~ TT^ tt'^) — m{K~TT^) is shown in Fig. ^a). A fit 
using a double Gaussian for the signal and the function 
shown in Eq. ^ for the background yields 69 000 ± 450 
signal events (statistical error only). 

The candidates are then combined with all pho- 
ton candidates producing the distribution of the mass 
difference Am(Z)"7) — m{K^TT^j) — m{K~Tr~^) shown 
in Fig. EJb). In this case, the peak corresponding to 
jj*o _5, jjo^ signal is close to a large bump arising from 
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TABLE II: Summary of the results. The first errors are sta- 
tistical; the second represent systematic uncertainties. 



r{D:+ 


^ z?+^»)/r(7?r 




0.062 ± 


0.005 ± 


0.006 


b{d:+ 






0.059 ± 


0.004 ± 


0.006 


B{D'+ 






0.942 ± 


0.004 ± 


0.006 


r(D*o 


^ D"7r")/r(D*° 




1.74 ± 


0.02 ± 


0.13 








0.635 ± 


0.003 ± 


0.017 








0.365 ± 


0.003 ± 


0.017 



the reflection of D*^ D'^tt^ in which one photon is pro- 
duced by tt'^ decay (the same reflection appears in D*^ 
decay but with a lower rate and less distinctive shape). 
Most of this bump is avoided by limiting the analysis to 
Am > 95 MeV/ . The remainder of the background is 
modeled using the function 



/3(Am) = Ar(^l + exp(^- 
X (Ato^ -f aAm 4 



Am — m(7r'^) 



(3) 



(Note that this function is similar to that of Eq. but 
differs in the sign of the exponential term.) The signal is 
modeled by the Crystal Ball function (Eq.|2Jl. The result- 
ing fitted signal consists of 67 880 ±670 events (statistical 
errors only). 

Efficiencies and systematic uncertainties are deter- 
mined using the procedures described for r(D*"'" 
D+TT°)/T{D*+ D+-i). Efficiencies of e{D°TT^) = 
0.037±0.002 and e{D^"f) = 0.064±0.002 and an efliciency 
ratio of e{D^TT^)/e{D^-i) ^ 0.58 ± 0.03 are found. The 
latter is consistent with the value of e(£'^7r°)/e(Z3+7). 
The ratio T{D*^ D°tt°)/T{D*° D°-f) = 1.74 ± 
0.02 (stat.) ± 0.13 (syst.) is obtained. 

The branching ratio measurements are summarized in 
Table Hll By assuming that the D*"*" meson decays only 
to DfiT^ and and that the D*° meson decays only 



to _D°7r° and D^j, it is possible to calculate the branching 
fractions, which are also listed in Table UTI 

In summary, the branching ratio r(I?*+ D^n^)/ 
T{D*+ D+7) = 0.062 ± 0.005 (stat.) ± 0.006 (syst.) 
has been measured and is consistent with the previous 
measurement 0, but has higher precision. Also deter- 
mined is the ratio T{D*^ -> D"Tr°)/T{D*" D^-f) = 
1.74±0.02(stat.)±0.13(syst.). This resuh is in agreement 
with, but is more precise than, the world average 

It has been proposed that the decay D*"*" Dfn° pro- 
ceeds via 77 — 7r° mixing and calculations based on Chiral 
perturbation theory Q predict B{D*'^ Dfir'^) « 1- 
3% based on current measurements of B{D*^ D'^j) = 
1.6 ± 0.4% |8j. Newer theoretical estimates in a relativis- 
tic quark model predict B{D*+ -> D+tt°) w 13%, 
somewhat larger than our measurement. 
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